arXiv:1504.00684vl [astro-ph.HE] 2 Apr 2015 


Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 6 April 2015 (MN DTpX style file v2.2) 


Study of the reflection spectrum of the accreting neutron star 
GX 3+1 using XMM-Newton and INTEGRAL 


Pintore F. \ Di Salvo T. 2 , Bozzo E. 3 , Sanna A. 1 , Burderi L. 1 , D’Ai A. 2 , Riggio A. 1 , Scarano F. 1 , Iaria R . 2 

1 Universita degli Studi di Cagliari, Dipartimento di Fisica, SP Monserrato-Sestn, KM 0. 7, 09042 Monserrato, Italy 

2 Dipartimento di Fisica e Chimica, Universita di Palermo, via Archirafi 36 - 90123 Palermo, Italy 

3 ISDC - Science data center for Astrophysics, Ch. d’Ecogia 16, 1290, Versoix (Geneva), Switzerland 


6 April 2015 


ABSTRACT 

Broad emission features of abundant chemical elements, such as Iron, are commonly seen in 
the X-ray spectra of accreting compact objects and their studies can provide useful informa¬ 
tion about the geometry of the accretion processes. In this work, we focus our attention on 
GX 3+1, a bright, persistent accreting low mass X-ray binary, classified as an atoll source. 
Its spectrum is well described by an accretion disc plus a stable comptonizing, optically thick 
corona which dominates the X-ray emission in the 0.3-20 keV energy band. In addition, four 
broad emission lines are found and we associate them with reflection of hard photons from 
the inner regions of the accretion disc where doppler and relativistic effects are important. We 
used self-consistent reflection models to fit the spectra of the 2010 XMM-Newton observation 
and the stacking of the whole datasets of 2010 INTEGRAL observations. We conclude that the 
spectra are consistent with reflection produced at ~ 10 gravitational radii by an accretion disc 
with an ionization parameter of £ ~ 600 erg cm s -1 and viewed under an inclination angle of 
the system of ~ 35°. Furthermore, we detected for the first time for GX 3+1, the presence of 
a powerlaw component dominant at energies higher than 20 keV, possibly associated with an 
optically thin component of non-thermal electrons. 

Key words: accretion, accretion discs - X-rays: binaries - X-Rays: galaxies - X-rays: indi¬ 
viduals 


1 INTRODUCTION 


Broad, and possibly asymmetrical, emission features are commonly 
fou nd in the X-r ay spe ctra of accreting Galactic black holes (BHs; 


Miller et at e.g. 2002!), supermassive BHs in Ac tive Galactic Nu¬ 
clei ( AGN: lTanaka~a l . 1995; FabianetalJ|2002 ) and neutron star 
systems INS: e.g. |Bhattacharwa & Strohmaveil 20071 : iD’Ai et aD 
120091 : IdTsalvo et alj|2009h . The main observed reflection features 
can be ascribed to transitions of the K-shell of iron in the en¬ 
ergy range 6.4-7.0 keV, according to its different state of ioniza¬ 
tions. In addition, broad features associated to lighter chemical 
elements have b een found in the sp e ctra of a number of accret¬ 
ing sources (e.g. [di Salvo et al JI20091 : llaria et afll2009l : iD’Ar et aD 


120101 : lEgron et al. 2013). It has been proposed that the broaden¬ 


ing of such features can be produced by either i) relativistic reflec- 
tion of hard photons from the surface of the accretion disc (e.g . 
Fabianet aD 1 19891: iRevnolds & Nowakl 120031: Fabian et al] 1200*1 


Mattll2006h or ii) Compton broadening if the line is emitted in the 


inner parts of a moderately, optically thick accretion disc corona 
(ADC), created by the evaporation of oute r layers of the dis c 
[White & Hold! 19821: iKallman & White!! 19891: IVrtilek et al.ll 19931) . 
iii) Another possibility is Compton down-scattering due to a narrow 
wind shell ejected at mildly relativistic velocities at some disc radii 


where the local radia tion force overco mes the local d isc gravity 
[Titarchuk et al.l20Qp) . We note also that lNg et~afl [20 1 (1) proposed 
that in many cases the asymmetrical shape of the spectral features 
seen in XMM-NewtonfEPlC-pn data may be due t o a wron g esti¬ 
mation of pile-up, although, on the contrary, iMiller et~aD 120101) 
showed that modest level of pile-up should not significantly affect 
the spectral analysis. However, the relativistic reflection scenario 
described above appears the more suitable explanation in many 
cases (e.g. Di Salvo et al. 2015, in press). 

The broad lines produced by relativistic reflection of hard pho¬ 
tons from the disc are usually double-peaked due to the motion of 
matter of the disc. Furthermore, the blueshifted peak is also rela- 
tivistically boosted as a consequence of the matter in the accretion 
disc approaching the observer at speed values of the order a few 
tenths of the speed of light when it is at a few gravitational radii 
from the compact object. The broad lines are then gravitationally 
redshifted when produced close to the compact object. The combi- 
nati on of s uch effects highly distorts the spectral shape of the line 
[Fabian et alj 19891) . These effects are observed in many broad lines 
of abundant chemical elements as Iron. In addition to the broad Iron 
emission feature, commonly detected at 6.4—7.0 keV, an absorption 
edge at ~ 8 — 9 keV is usually observed. 
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Moreover, a large fraction of hard X-ray photons is reflected 
by the accretion disc producing a reflection continuum, due to 
the direct Compton scattering, commonly known as “Com pton 
bump” which peaks at around 20-30 keV (e.g. iGeorge & Fabianl 
1 199 ill . Even though often described independently, the discrete 
emission features and the reflection continuum are closely re¬ 
lated. Self-consistent models can instead investigate the physics 
of the reflection processes takin g into account both spectral com- 
ponents at the same time jjRoss^t dll 19991: iRoss & FabiaiJl2007l: 


Reisetah||2008k Cackett et alJbOloT D’Ai et al.ll20ld : Egron et al. 

20131 : Pi Salvo et alJ 2015l) ~The detailed study of the reflection 


component plays therefore a fundamental role in the investigation 
of the properties of the compact objects (BHs and NSs) and the 
physical mechanisms of the accreting material around them. 

In this work we investigate the spectra l properties of the 
Galactic source GX 3+1 IBowver et al .1119651) . which is a bright, 
persistent accreting low mass x-ray binary (LMXB). A number 
of type I X-ray bursts, likely due to burning helium, has been 
detected, with de cay timescales ranging from tens of seconds 
up to hours (e.g. Makishi ma et al. 1983; jdenHartog_et_aL 20031 


Kuulkers & van der Klisl l2000l: iKuulkersI 120021: IChenevez et al] 


2006h . The presence of thermonuclear bursts suggest that the com¬ 


pact object hosted in GX 3+1 is an accreting NS. The_estimated 
upper limit to its distance is ~ 6.1 kpc jKuulkers & van der Klisl 
l2000h . A probable quasi periodic oscillation (QPO) with a fre¬ 
quency of 8 Hz has been detected during an EXOSAT observation 
iLewin et all 19871) . Based on its variabilit y and spectral properties, 
GX 3 + 1 was classified as an atoll source IHasinger & van der Klisl 
Il989h . Ivan den Berg et al] d2014l) identified a counterpart of GX 
3+1 in the near-infrared (NIR) with a Ks=15.8 +/- 0.1 mag suggest¬ 
ing that the donor star is not a late-type giant and that the NIR emis¬ 
sion is dominated by the heated outer accretion disc. Spectral anal¬ 
yses of the source showed that its X-ray spectrum can be described 
by a model comprising of a soft blackbody component, most likely 
associated with the accretion disc, and a comptonizing component, 
produced by an optically thick electron population located close to 
the NS boundary layer 1 Mainardi et ali 2010|^ Piraino et~ai1 1 201 2l ; 


ISeifina & Titarchukl l2012h r Seifina & Titarchuk d2012h T~ using all 

the Rossi X-ray Timing Explorer (RXTE) and BeppoSAX data, 
showed that the comptonizing component is quite stable suggesting 
that the energy release in the transition layer between the accretion 
disc and the NS dominates over the energy release in the accre- 
tio n disc. A number of broad emission lines were also reported by 
Ipiraino et alJ d2012h and associated with the Ka transitions of Ar, 
Ca and Fe. These features have been ascribed to relativistic reflec¬ 
tion and the inferred inner disc radi us of ~3 0 gravitational radii 
indicated a truncated accretion disc JPiraino et all 201 3) . 

In next section, we will carry out a deep spectral analysis 
of the broadband spectrum of the stacking of the whole 2010 IN¬ 
TEGRAL observations and the 2010 XMM-Newton observation. In 
particular, we will study the prop erties of the reflection component, 
as shown in lPiraino et al.l d2012f) . and extending it with the use of 
self-consistent models. The paper is structured as follows: in Sec¬ 
tion [2] we summarize the data selection and reduction procedures, 
in Section[3]we present the spectral results of GX 3+1 and discuss 
them in Section[4] 


2 DATA REDUCTION 

The datasets considered here comprised the XMM-Newton obser¬ 
vation (Obs.ID. 0655330201), taken in Timing mode on September 


1st, 2010, combined with the whole available 2010 INTEGRAL ob¬ 
servations. 

XMM-NewtonfEPIC-pn data were reduced using the calibra¬ 
tion files (at the date of April 25th, 2014) and the Science Analysis 
Software (SAS) v. 13.5.0. At high count rates. X-ray loading and 
charge transfer inefficiency (CTI) effects have to be considered be¬ 
cause they modify the spectral shape and produce an energy shift 
on the spectral features (e.g. Guainazzi et al. 2008, XMM-CAL- 
SRN-24f0). We corrected these effects by adopting the rate de¬ 
pendent pulse height amplitude corrections (RDPHA) which have 
been shown to pro vide a better accuracy of the energy scale calibra¬ 
tion JPintore et alj2014h . We extracted the spectra from events with 
PATTERN^ 4 (which allows for single and double pixel events) and 
we set ‘FLAG=0’, retaining events optimally calibrated for spectral 
analysis. Source and background spectra were then extracted se¬ 
lecting the ranges RAWX=[31:43] and RAWX=[3:5], respectively. 
EPIC-pn spectra were rebinned in energy with an oversample of 
7 using the specgroup task. We looked for pile-up in the EPIC-pn 
data by comparing the source spectrum obtained considering all the 
columns of the CCD in the range RAWX=[31:43] and the spectra 
obtained excising one, three and five central brightest columns of 
the aforementioned RAWX range. We found that pile-up was only 
marginally important and it could be clearly mitigated removing 
only the brightest central column of the EPIC-pn CCD. Hereafter 
we will refer, if not specified, to spectral analysis without the cen¬ 
tral column. 

We extracted first and second order RGS spectra using the 
standard rgsproc task, filtered for periods of high background. After 
checking for spectral consistency of RGS 1 and 2 in both first and 
second order, we stacked together the RGS 1 and RGS2 data adopt¬ 
ing the tool rgscombine which combines, separately, first and sec¬ 
ond order spectra. We did not find any feature in the RGS datasets, 
therefore we grouped the two final spectra with a minimum of 100 
counts per energy channel. Since the background, in both orders, 
was not negligible for energies lower than 0.9 keV, we analyzed the 
RGS spectra in the range 0.9-2.0 keV. Even though the mean count 
rate is at least a factor of 3 lower (RGSl ~3.4 and RGS2 ~ 1.7 cs 
s -1 ) than the nominal RGS pile-up threshold, we investigated our 
data looking up for pile-up effects and we found that pile-up does 
not significantly affect our data because we verified that first and 
second order spectra are in complete agreement with each other. 
This allows us to confirm that pile-up problems are negligible for 
this source. 

For the scope of the present paper, we also considered all 
the available INTEGRAL data collected in 2010 in the direction 
of GX3+1. The dataset included “science windows” (SCWs), i.e. 
the INTEGRAL pointings with typical durations of ~2 ks each, 
from satellite revolution 895 to 981. All data were analyzed by 
using vers ion 10.0 of the OS A software distributed by the ISDC 
dCourvoisier et alj |2003h. We used dat a from IBIS/ISGRI (17-80 

ir 


keV, iLebrun et al.ll2003l : lubertini et al.ll2003h in which the source 
was located to within 12° from the center of the instrument field of 
view (FoV), thus avoiding problems with calibration uncertainties. 
We also extracted JEM-X1 and JEM-X2 ( Lund et ah] 120031) data 
from the same set of SCWs. The ISGRI (JEM-X) response matri¬ 
ces were rebinned in order to have 37 (32) energy bins spanning the 
range 20-180 keV (3-35 keV) for all spectra. This optimized the 
signal-to-noise (S/N) of the data. 

We fitted simultaneously the XMM-Newton (RGS and EPIC- 


http://xmm2.esac.esa.int/docs/documents/CAL-SRN-0248-l-0.ps.gz 
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pn) and INTEGRAL (JEM-X1, JEM-X2 and INTEGRAL/ISGRI) 
spectra, using XSPEC V. 12.8.2 lAmaudlll996h . We selected the en¬ 
ergy ranges 0.9-2.0 keV, 2.4-10.0 keV, 10.0-20.0 keV and 20.0- 
50.0 keV, for RGS, EPIC-pn, JEM-X1, JEM-X2 and ISGRI, re¬ 
spectively. We ignored the EPIC-pn channels at energies lower than 
2.4 keV as we found a soft excess which was not reconciled with 
the RGS data, suggesting that the EPIC-pn calibrations are still un¬ 
certain below this ener gy (i nternal XMM-Newton report CAL-TN- 
008:0; see also lPiraino et al.ll2012i ). 

We introduced a multiplicative constant model in each fit in 
order to take into account the diverse flux calibrations of the con¬ 
sidered instruments and possible differences on the flux caused by 
non simultaneity of XMM-Newton and INTEGRAL observations. 
We fixed to 1 the constant for the EPIC-pn spectrum and allowed 
the other constants to vary. In general, the range of this parame¬ 
ter is consistent within 10% for the RGS instruments and does not 
vary more than ~ 20% for the INTEGRAL/ISGRI and IEM-X, in 
comparison with the EPIC-pn constant. 


3 SPECTRAL ANALYSIS 

The aim of the spectral analysis presented here is to provide a char¬ 
acterization of the X-ray emission from GX 3+1 in a broadband 
energy range. We are particularly interested in the reflection com¬ 
ponent whose properties can be accurately described only when 
the overall continuum is well constrained. Previous works (e.g. 
Seifin a~& Titarchukll20l2h have shown that the high energy com¬ 
ponent of GX 3+1 is weakly variable. 

We analyse the 2010 XMM-Newton observation with a net ex¬ 
posure time of 55 ks. In Figure [I] we show the total light curve in 
the 0.3-10 keV energy band: during the exposure, a type I burst, 
which characterizes the temporal properties of the source, was ob¬ 
served and GX 3+1 presented also an unambiguous evidence of 
flux variability on timescales of ~ 5 ks. To test the contribution 
of the burst in the spectral parameters, we fit RGS and EPIC-pn 
spectra (with and without the burst time interval) adopting a self- 
consistent model, as showed in section [3711 and we find that the 
X-ray burst does not affect our spectral parameters. For the tempo¬ 
ral variability, we instead investigate if such variability could also 
drive spectral changes. To test this, we follow two alternative ap¬ 
proaches based both on temporal and flux variability: in the first 
one, we split the observation in 11 temporal segments of 5 ks each; 
while, in the second one, in order to further improve the signal- 
to-noise ratio, we divide the observation in three flux binned seg¬ 
ments (low flux (< 370 cs s -1 ), medium flux (between 370 cs s -1 
and 400 cs s -1 ) and high flux (> 400 cs s -1 )) and with a similar 
amount of collected counts. We extract an EPIC-pn spectrum, pile- 
up corrected, for each segment and we fit the obtained spectra in 
the 2.4-10 keV energy range, where the calibration of the instru¬ 
ment is more accurate. We ad opt a simple comptonization model 
based on an NTHCOMP jZdziarski et alJl99fih plus a relativistic line 
(DISKLINE. lFabian et al.l 19891 ; see below for more details). We find 
that all the spectral parameters inferred by the spectral analysis are 
consistent within the uncertainties (at 3 a confidence level). We thus 
report in the following only on the average spectrum of the whole 
XMM-Newton observation. 

We note that the Monitor of All-sky X-ray Image (MAXI) 
telescope observed GX 3+1 in a quiete stationary flux level in 

2 http://xmm2.esac.esa.int/docs/documents/CAL-TN-0083.pdf 


GX 3+1 



Figure 1. XMM-Newton! EPIC-pn lightcurve in the 0.3-10 keV energy 
range, corrected for pile-up excising the brightest central column and sam¬ 
pled with bins of 20 seconds. A type-I burst and clear time variability are 
observed. 
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Figure 2. Unfolded E 2 f(E) XMM-NewtonlEPIC-pn (black), XMM- 
Newtonl RGS spectra ( red and green), INTEGRAL! JEM-X1 and X2 (cyan 
and magenta), and INTEGRAL/ISGRI (blue). The dashed curves represent 
the BBODY and the NTHCOMP components (see text). For display purpose 
only, RGS data have been rebinned with a minimum of 20rr and maximum 
number of energy channels of 20. 

the 10-20 keV energy band during the entire 2010, with average 
variability lower than a factor of 2. Therefore, we also decide to 
stack all the available 2010 INTEGRAL observations of GX 3+1 
in order to improve the signal-to-noise ratio in the band 10-50 
keV. In addition, this choice allows us to collect very high quality 
data for a period which is reasonably close to the XMM-Newton 
observation. However, although the flux variability was limited 
during the 2010, we verify that also the spectral variability was 
negligible during the year: we compare the INTEGRAL spectra 
(both ISGRI and JEM-X instruments) obtained by stacking the 
data of the whole year and those obtained by stacking the data 
of about 100 days centered at the XMM-Newton observation. We 
select 100 days as a test in order to have a good quality counting 
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Table 1 . Best fit spectral parameters obtained with the absorbed continuum BBODY+NTHCOMP+POWERLAW model. Column 1, 2, 3 and 4 refer to best fit 
models where the reflection component is described by GAUSSIAN, DISKLINE, RFXCONV or REFLIONX models, respectively. Errors are at 90% for each 
parameter. 


Model 

Component 

( 1 ) 

( 2 ) 

(3) 

(4) 

PHABS 

Nj, (10 22 cm- 2 )“ 

2 19+ 0 - 08 
z,iy -0.09 

2 l + 01 

2 oi + 0 - 03 

z - ul -0.04 

9 Q9+0.03 
z,uz -0 04 

BBODY 

kT b( , (keV ) 6 

0-54 ±g;Sl 

0 57 + 0 - 02 
u, °'-0.03 

U.OU _ 0 03 



Norm. <10 3 ) c 

20 - 2 ^o 2 

17 q + 0 - 2 

11 ' y —0.3 

c c+0.9 

0.9 

10 . 2±°4 

NTHCOMP 

kT e (keV) d 

0 q+ 0.6 

O-U -0.3 

q 1 + 0.8 

^• a -0.5 

9 oe :+ 0 - 04 
z,zo -0.04 

9 9-1 +0.04 
z - zl -0.04 


r e 

9 y+ 0.2 
^•'-0.3 

9 9 + 0.3 
Z,y -0.5 

i ' 79 +0.03 
- L - < 1 - 0.02 

l. 8 U_o 05 


kV 

1 9+0.07 
i - z -0.07 

1 26 +0 ' 07 

0.08 

0.69l“;“ 6 

0 92 +0 ' 10 

0.08 


Norm. <10 2 )9 


io±i 

91+ 2 

Z 1_2 

12“*" 2 

iz _ 2 

POWERLAW 

r e 

4 -+° 0 .i 

4 o + 0 - 6 
^• u -0.7 

3 36+ 010 

O.OO-0 05 

q q + 0.3 

°-°-o.i 


Norm. 9 

2 I + 0,4 
Z,i - 0.2 

1 5+ 0 - 7 
i, 0 -0.4 

1 2" 1 " 0 ' 1 
i, ^- 0.2 

u,y - 0.1 

GAUSSIAN 

Energy (keV) h 

9 zq+ 0-03 
^•Oo _ 0 Q 2 

_ 

- 

_ 


G i 

0 06 + ° 05 
u.uo _ 0 04 

- 

- 

- 


Norm. (10“ 3 y 

0 7 + 0,5 
u -'- 0.2 

- 

- 

- 

GAUSSIAN 

Energy (keV)^ 

q 09 +O.O 2 
' — 0.02 

_ 

- 

_ 


(7 ® 

0 m+ 0.03 
u * A ^-0.03 

- 

- 

- 


Norm. (10 _3 y 

1 3 +0 - 4 

1 O -0.3 


- 

- 

GAUSSIAN 

Energy (keV)^ 

0 qo+ 0.02 
J ' yo - 0.02 

_ 

- 

_ 


G i 


- 

- 

- 


Norm. (10 _3 y 

0.8+°;4 

- 

- 

- 

GAUSSIAN 

Energy (keV) h 

(■ /ro+ 0.02 

b-o 8_ 0 02 

_ 

- 

_ 


G * 

0 32 +0-02 
' j ‘ ja —0.02 

- 

- 

- 


Norm. (10 _3 P 

7 4+ 0 - 2 
^— 0.2 

" 

- 

- 

DISKLINE 

Energy (keV) h 

_ 

2 64 +0 ' 02 

0.02 

_ 

_ 


Norm. (10 _3 y 

- 

1 6 +0 - 4 

A '°—0.4 

- 

- 

DISKLINE 

Energy (keV) h 

_ 

q q 9 + 0.02 
J - J ^-0.02 

0 09 + 0.02 
J ' JZ ‘—0.02 

q qc+ 0.02 

0.02 


Norm. (10 _3 P 


1 5 "*" 0,2 

1 ■ 0.2 

1 9 +O .2 

1 ,z '— 0.2 

0 - 8 ±g? 

DISKLINE 

Energy (keV) h 

- 

3 qc+ 0.02 
- 0.02 

0 qo+ 0.02 

0.02 

q Q7+0-0 2 
' — 0.02 


Norm. (10 -3 y 


1 i+ 0.2 

1 ' 1 — 0.2 

Q q + 0.2 
uy -0.2 

0 6 +01 
u '°- 0.1 

DISKLINE 

Energy (keV) h 

- 

f. ^c+ 0.02 

6.65_o 02 

- 

- 


Norm. (10 _3 y 


9 c+ 0.2 
^ ,J - 0.1 

~ 

- 

RDBLUR 

BetorlO k 

- 

-2 63" 1 " 0 ' 05 

z.oj _ 0 05 

9 qn+0.06 

-Z.jU— o.07 



Rin (R< 7 / 

- 

32t;° 




R out (Rg) 771 

- 

10 4 

10 4 

10 4 


Inclination (degree) 71 


41+ 4 

41 —3 

35 4+ 0 - 7 


RFXCONV 

Fe/solar° 

_ 

_ 

=1 (fixed) 

_ 


Log^P 

- 

- 

9 ' 79 + 0.05 

Z -‘ ‘ -0.05 

- 


Rel refl (sr) 9 

- 


n 99 +O.O 5 
u ZZ -0.03 


HIGHECUT 

CutoffE 7- 

- 

- 

- 

=0.2 (fixed) 


FoldE 5 

- 

- 

- 

=2.7 x kT e 

REFLIONX 

Fe/solar° 

- 

- 

- 

=1 (fixed) 


Log^ 

- 

- 

- 

2 94 + 0 - 09 


Norm. (10 -6 ) 1 

- 

- 

- 

19 9 + 0.5 
lz - z -0.4 


x 2 /dof u 

2126.22/1674 

2091.52/1675 

2132.25/1677 

2142.65/1677 


a Neutral column density; b Blackbody temperature; c Normalization of the BBODY in unity of L 39 / D^ okpc , where L 39 is the luminosity in unity of 10 39 
erg s —1 and D\Qk pc is the distance in unity of 10 kpc; d Electrons temperature of the corona; e Photon index; f Seed photons temperature; 9 Normalization 
of NTHCOMP; h Energy of the line; 1 Broadening of the line; 9 Normalization of the line; k Power law dependence of emissivity; 1 Inner radius in units of 
gravitational radii ; 771 Outer radius in units of gravitational radii R g ; 71 Inclination angle of the binary system; 0 Ratio of Iron and Hydrogen abundance; p 
Logarithm of the ionization parameter, with the latter expressed in units of erg cm s —1 ; q Solid angle (0/27r) subtended by the reflector as seen from the 
corona; r high energy cut-off of the reflection component, set as 2.7 times the electron temperature of the comptonisation model; s folding energy of the 
cut-off; 1 Normalization of the REFLIONX component; u reduced x 2 ; 
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statistics. No significant spectral variability is detected. The 
total exposure times of the INTEGRAL instruments is then 700 
ks, 64 ks and 390 ks for ISGRI, JEM-X1 and JEM-X2, respectively. 

We initially investigate the spectral properties of GX 3+1 
adopting a model consisting of an absorbed (PHABS in XSPEC) 
blackbody (BBODY) plus a comptonizing component (NTHCOMP), 
as it was shown in previous works to b e a good representa¬ 
tion of the source X-ray continuum fe.g. | Main ardi et al.j 1 20 id : 
ISeifina & Titarchukl 120121 : IPiraino et al JI 201 We find a neutral 
col umn density (Ng) of ~ 1.6 x 10 22 cm -2 , using the abundances 
of lAnders & Grevessel il989l) . At soft energies, a blackbody 
component with a temperature (kTj,b) of ~ 0.6 keV is observed 
and it can be associated to a hot accretion disc as the inferred 
emission radius is 28 ± 2 km, which allows us to exclude emission 
from the surface of a 1.4 Mq NS. The comptonizing component 
can be instead described by an electron temperature (kT e ) of ~ 2 
keV, a seed photons temperature (kTo) of ~ 0.2 keV and a photon 
index of 1.4. This corona probably originates in the inner regions 
of the disc or more likely in the boundary layer. This model 
provides a poor fit (x 2 /do/=4881.49/1688; Fig. 0, although 
we note that it should be the best choice for the continuum. 
Indeed other similar spectral continuum models, for example 
as BLACKBODY+POWERLAW or a cold DISKBB l lMitsuda etaD 
119841) plus a hot BLACKBODY, were tried but they gave sig¬ 
nificantly worse statistical results (x 2 /dof= 112812.0/1700 and 
X 2 /dof= 5063.76/1690, respectively). However, the poor fit of 
the BB+NTHCOMP model is not surprising because residuals 
from the fit are observed at 2.6, 3.3, 3.9 and 6.4 keV, due to 
not modelled emission lines. We also notice the presence of an 
excess at energies higher than 15 keV. In order to improve the 
fit, we firstly introduce four GAUSSIAN components letting the 
broadening (cr) to vary freely amongst the lines. An improvement 
in the fit is found (x 2 /dof= 2261/1676) and the measured centroid 
energy lines is 2.71, 3.33, 3.97 and 6.66 keV. These lines can 
thus be associated with S XVI, Ar XVIII, Ca XX (or XIX) and 
Fe XXV. The emission features are broad, with widths comprise 
between ~0.1 and 0.3 keV (Fig. [3] top-left). This confirms the 
results of lPiraino et alJ l2012i ). although we note that we have been 
able to identify also the S XVI line, which is only marginally seen 
in the spectral residuals of their work but statistically significant 
in our (Ay 2 =71 for 2 additional dof). Furthermore, we add that 
the NTHCOMP parameters are now changed: the photon index, 
the electron temperature and the seed photons temperature are 
increased to 2.7, 3 keV and 1.2 keV, respectively (see TableQJ. 


The presence of an excess in the residuals from the fit at en¬ 
ergies lower than 15 keV was not affected by the introduction of 
the lines. Therefore this may be the evidence that it is produced 
by an additional spectral component that we tentatively model 
with a PO WERLAW. Hard powerlaws are not unusual in accret- 
ing NS (e.g Dp Salvo et al .1200(1 1200 1 1 : ITari a et al.l 200 ll llaria et al l 


200_Jy_ di_Salvo_t+af 20041 : Paizis et al. 20061 : Piraino et al. 1 120071; 
Migli ari et al.ll2007l) . although it has never been clearly observed 
in the spectra of GX 3+1. The introduction of this component sig¬ 
nificantly improves the fit (fl 2 / dof=2126. 22/1674) and, hereafter, 
it will be considered in all the spectral fits as its presence is al¬ 
ways statistically significant (Ax ' 2 > 70 for 2 additional degree of 
freedom). We find that the photon index is very steep (4.7 ± 0.5) 
and it is not able to completely flatten the high energy excess. We 
note however that a fit, with a slightly worse statistical significance 


(Ay 2 


10 ), can also be found, when the photon index is con¬ 


sistent with ~ 3.4. Although the former is statistically better, the 
photon index of the latter may be more likely. 

We subsequently substitute the four GAUSSIAN models with 
relativistic broadened emission features (DISKLINE in XSPEC), in 
order to test the goodness of the relativistic reflection interpretation. 
For all the lines, we kept tied the inner disc radius (R;„), the outer 
disc radius (R out), the inclination angle and the emissivity index of 
the illuminating source (betor 10). We remark that the spectral fit is 
insensitive to the outer disc radius parameter, hence we left it frozen 
to the value of 10 4 gravitational radii (R 9 = 2GA/*/c 2 , where 
M* is the NS mass). These models provide a significant improve¬ 
ment if compared with the simple GAUSSIAN models (x 2 /dof = 
2091.52/1675) as would be expected if the lines are produced as a 
result of reflection (see Section |T}. We found that the energy of the 
lines are consistent with those previously determined with the gaus- 
sian models and the relativistic parameters indicate that the system 
is inclined with an angle of ~ 41°. In addition, our results sug¬ 
gest that the inner accretion disc is possibly truncated away from 
the innermost^stable circular orbit (R,„ ~ 31 R g ), confirming the 
findings of IPiraino et al.1 d2012l) . However, we propose that for a 
truncation radius of ~30 R g , we would expect an inner disc tem¬ 
perature of ~ 0.3 keV which significantly differs from the 0.5 keV 
temperature found by the spectral analysis and, in particular, to its 
corresponding emission radius of 22 ± 2 km, (assuming a distance 
of 6.1 kpc), i.e. ~11 R 9 . We note how this result may be affected by 
the simplicity of the spectral model and it suggests that the adopted 
reflection model cannot provide a coherent picture for the under¬ 
lying physics of the accretion processes in GX 3+1. Therefore a 
self-consistent model appears to be the best choice for the spectral 
analysis of the source. 


3.1 Self-consistent models 

The existence of four relativistic emission features suggests that 
they are produced by reflection of hard photons from the sur¬ 
face of the accretion disc. This makes GX 3+1 one of the most 
intriguing sources to be studied adopting self-consistent models 
which describe the reflection continuum component and the dis¬ 
crete features in accordance with the chemical composition of 
the accretion disc. To accomplish this study, we make use of a 
commo nly used self-consistent m odel: the RFXCONV convolution 
model I Kolehma inen et aPl201 ll) . RFXCONV depends on five pa¬ 
rameters: the relative reflection normalization rel re /;, the redshift 
of the source, the iron abundance relative to Solar (Fe/Solar), 
the cosine of the inclination angle and the ionization parameter 
(Logtf) of the disc reflecting surface. We convolve this model with 
the NTHCOMP and the POWERLAW models, as they are the contin¬ 
uum source of hard photons of the spectrum. We further convolve 
the RFXCONV component with a relativistic kernel (RDBLUR in 
XSPEC) to take into account relativistic distortion due to a rotating 
disc close to a compact object. As RFXCONV does not take into ac¬ 
count the presence of Ar XVIII and Ca XIX features but only the Fe 
and S lines (which are more abundant), we include two DISKLINE 
components for them. In particular, we link the parameters of RD¬ 
BLUR to those of the DISKLINE. We note that the Iron line is the 
strongest emission feature in the spectra and it drives the spectral 
parameters of the other lines. However, when the lines are left free 
to vary independently, the spectral parameters converge towards 
values consistent within the errors. It suggests that all the lines have 
a common origin. The best fit with this model does not introduce 
any significant statistical improvement (f 2 /dof =2132.25/1671, 
see Table QJ-column #3 and figure \3f bottom-left) with respect to 
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phabs*(bbody+gaussian+gaussian+gaussian+gaussian+nthcomp+powerlaw) 


phabs*(bbody+diskline+diskline+diskline+diskline+nthcomp+powerlaw) 



12 5 10 20 

Energy (keV) 

phabs*(bbody+diskline+diskline+rdblur*rfxconv*(nthcomp+powerlaw)) 



12 5 10 20 

Energy (keV) 

phabs*(bbody+diskline+diskline+highecut*rdblur*reflionx+nthcomp+powerlaw) 



Figure 3. Unfolded E 2 f(E) XMM-Newton/EPIC-pn {black), XMM-Newton/RGS spectra {red, first order, and green, secondo order), INTEGRAL! JEM-X1 
and X2 {cyan and magenta) and I NTF, G Aft / VIS G R i {blue). The top panel shows the best fit with the GAUSSIAN (top-left) or DISKLINE (top-right) models; the 
bottom panel shows the best fit with the RFXCONV (bottom-left) or REFLIONX (bottom-right) self-consistent models. The dashed curves represent the BBODY 
{light blue), the NTHCOMP (light red), the POWERLAW (dark grey), the GAUSSIAN or DISKLINE (purple ) and the self-consistent (orange) components. For 
display purpose only, RGS data have been rebinned with a minimum of 20a- and maximum number of energy channels of 20. 


the DISKLINE model but allows us to give a more physical inter¬ 
pretation of the spectrum of GX 3+1. We find that the solid an¬ 
gle (tl/2n) subtended by the reflector as seen from the c orona is 
0.22 ± 0.5 steradians, confirming the expectations of iDove et al .1 
Jl997t) . which estimated a solid angle of 0.3 for a spherical geom¬ 
etry. We find also that the logarithm of the ionization parameter of 
the disc is ~ 2.8, which can widely explain the ionization state of 
the four emission lines observed in the spectrum. The inclination 
angle of the system is found to be consistent with 35°. The iron 
abundance is compatible with the solar value (assumed to be 1) as 
the fit was almost insensitive to this parameter and, for this rea¬ 
son, we kept it fixed to unity. However the broadband continuum, 
with the introduction of this self-consistent model, is significantly 
changed: the photon index of the NTHCOMP is decreased to ~ 1.8 
as well as the electron temperature which is ~ 2.3 keV. In addition, 
the most significant change is seen in the POWERLAW photon index 
and seed photons temperature which are decreased towards ~ 3.4 
and ~ 0.7 keV, respectively. The photon index of the POWERLAW 


is now less extreme than that found in the previous section (> 4) 
and it is able to fit satisfactorily the high energy excess. This re¬ 
sult suggests that the continuum reflection component is required 
to model correctly the broadband continuum. In addition, the in¬ 
ferred inner disc radius is consistent with 7 R 9 and hence largely 
discrepant in comparison to the simple DISKLINE model. This lat¬ 
ter result is definitely more consistent with the temperature of the 
disc (i.e. ~0.5-0.6 keV) and with its corresponding emission ra¬ 
dius of ~ 15 ± 2 km, i.e. R g . Furthermore, we note that the 
seed photons temperature and the blackbody temperature are very 
similar, suggesting that the inner region of the disc is the source of 
seed photons for comptonization and that, more likely, the optically 
thick corona is the boundary layer. 

To validate our interpretation, we test ano ther widely used 
self-consistent model: the REFLIONX model jRoss & Fabianl2005h . 
The powerlaw photon index of the REFLIONX component is tied to 
the photon index of the NTHCOMP model, which appears to be the 
dominant source of hard photons (see Figure [3). We then convolve 
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the REFLIONX model with a relativistic kernel (RDBLUR, and also 
with a HIGHECUT model in order to limit the spectral extension of 
the REFLIONX model and mimic the roll-over of the comptoniza- 
tion component. The folding energy is fixed at 0.2 keV while the 
high energy cut-off is linked to the electron temperature of the NTH- 
COMP model. The parameters of the RDBLUR are then linked to 
those of the DISKLINE models used for the Ar XVIII and Ca XIX 
lines. 

Not surprisingly, the best fit with this model is statistically 
consistent with the previous one (x 2 /dof=2142. 65/1677, see Ta¬ 
ble □ column #4 and figure [3] bottom-right ), as well as the spec¬ 
tral findings. We note that the reflection component is produced 
by chemical elements with a logarithm of the ionization of ~2.9, 
reasonably consistent (within the uncertainties) with the RFXCONV 
estimate. The inclination angle of the system is estimated again to 
be 35°. The relativistic parameters are also in line with the previous 
ones, and we find an inner disc radius roughly consistent with the 
innermost stable circular orbit. This confirms that the RFXCONV 
and REFLIONX models give results consistent with each other and 
are clearly reproducing the same type of reflection process and ge¬ 
ometry of the accretion flow, reinforcing our interpretations and 
conclusions. 

It is important to mention that the self-consistent reflection 
models that we used in this work do not take into account the Ar and 
Ca lines which appear to be important in the spectrum of GX 3+1. 
For this reason, in order to further investigate the nature of the Ar 
and Ca lines, we also tried the RELXILL model dGarcia et all2QI4l . 
i.e. the relativistic version of the XILLVER model); this is currently 
the only self-consistent reflection model including Ar and Ca tran¬ 
sitions and is, therefore, the best model to fit these low energy lines. 
We tested the latest public version of this model (v0.2g) where Ar 
and Ca transitions are included but it is not possible to indepen¬ 
dently vary their element abundances. The residuals, with respect 
to the best-fit, still point to an unaccounted excess at the energies 
of the resonant lines, suggesting a possible overabundance of these 
elements. We find, from the DISKLINE models, that the ratio of the 
equivalent widths of the Ar (Ca) line with respect to that of the Fe is 
~10 (13)%. Using the simple calculations showed in lDi Salvo et al.l 
d2015h . we might expect an overabundance of these elements of ~4 
times the solar abundance. We note only that, in the GX 3+1 spec¬ 
trum, they are highly statistically significant and, in addition, these 
lines are very often found in XMM- Newto n sp ectra (and perhaps 
also in Suzaku spectra, see lDi Salvo et alj2015l) of bright LMXBs. 
In conclusion, at this moment we cannot infer any unambiguous 
indication about the abundance of these elements, which need to be 
further investigated with suitable reflection models. 

Finally, we discuss the importance of the POWERLAW compo¬ 
nent: we verified that, without such a component, the fits with self- 
consistent models would provide unphysical results. In particular, 
the BLACKBODY temperature would have been decreased toward 
lower temperatures (~0.15 keV) although Rwas still close to 6 
R s . Such a low temperature could only be consistent with a highly 
truncated accretion disc, which does not seem the case for GX 
3+1. For completeness, we mention that we also tried an alternative 
model of the continuum in which the seed photons for comptoniza- 
tion were considered as originating from an accretion disc and with 
a blackbody component associated t o a hot b oundary layer/NS sur¬ 
face emission (“Western” model, e.g. lwhite et all 19881) . The statis¬ 
tics of this model is marginally worse (y 2 /dof= 2157.16/1677) 
than the continuum model adopted above, but the results regarding 
the reflection parameters and the significance of the hard powerlaw 
are totally in line with those previously found. The only difference 


is related to the seed photons temperature (~ 0.5 keV) and the 
blackbody temperature (~ 0.9 keV) which are, however, simply 
reversed if compared with the previous ones. 


4 DISCUSSION AND CONCLUSIONS 

We analysed a XMM-Newton observation, taken in timing mode, 
with the aim to investigate the reflection properties of the persis¬ 
tent accreting NS GX 3+1. In addition, we also made use of all the 
available INTEGRAL observations collected from January to De¬ 
cember, 2010. The whole INTEGRAL data were stacked together 
to improve the signal-to-noise ratio. We have been able to study the 
X-ray emission from the source in the energy range 0.9-50 keV. GX 
3+1 is a bright, accreting NS located in the direction of the Galactic 
center. We estimated an unabsorbed total flux of ~ 7 x 10 -9 erg 
cm -2 s -1 and a luminosity of ~ 3 x 10 3 ' erg s -1 (for a distance of 
6.1 kpc). This corresponds to ~20% of the Eddington luminosity. 
The source is characterised by a complex spectral shape with asym¬ 
metrically broadened emission features which were associated to 
reflection of hard photons from the surface of the accretion disc. 
They were identified with the transitions of Iron, Calcium, Argon 
and Sulphur in H-like ionization state. We described their prop¬ 
erties adopting either simple GAUSSIAN or DISKLINE models, or 
self-consistent models. Our results wi th the DISK LINE models are 
consistent with the previous work of lPiraino et alj I 2 OI 2 I) . although 
the disc temperature is significantly lower in our work (~0.6 keV 
versus 0.9 keV) possibly due to the use of the RGS data. How¬ 
ever, such models do not provide a global picture for the accretion 
mechanisms on GX 3+1. On the other hand, the spectral param¬ 
eters obtained adopting self-consistent models are more coherent 
with each other and therefore, hereafter, we discuss the results of 
the self-consistent models as they provide a more physical interpre¬ 
tation of the X-ray emission from GX 3+1. 

We modeled the continuum with a blackbody, more likely as¬ 
sociated with the accretion disc, and a comptonizing component. 
The latter component is described by a cold (~ 2 keV) electron 
population having a relatively flat spectrum (r ~ 2) which de¬ 
creases rapidly above 20 keV. This component is dominant at en¬ 
ergies between 1.5-30 keV, corresponding to about 70% of the to¬ 
tal flux emission, and the stabilit y of this high energy comp onent 
has been widely investigated (e.g. lSeifina & Titarchukl(20121 . This 
makes feasible our choice to stack all the available INTEGRAL ob¬ 
servations. Such a component is probably produced in the transition 
layer around the NS and can par tially cover the inner regions of the 
disc iSeifina & Titarchukll2012i) . The accretion disc is instead rel¬ 
atively hot (~ 0.5-0.6 keV) and we inferred an emission radius of 
15 ± 2 km, suggesting it is close to the compact object, and with 
a temperature which is very similar to the seed photons tempera¬ 
ture of the optically thick corona, strongly supporting the location 
of the corona close to the NS. In addition, we identified a steep 
powerlaw component (T ~ 3.4) that does not show any indication 
of a high energy cut-off and that is predominant at energies higher 
than 30 keV. It is the first time that this component is observed in 
GX 3+1 and we found that it contributes to around 20% of the to¬ 
tal flux. Hard powerlaws, also known as hard tails, are commonly 
seen during both low/hard and high/soft states of accreting black 


holes (McClintock & Remillarc 

2006) and also in a number of ac- 

creting NS (e.g. 

Di Salvo et alj 

2000, 2001: D’Amico et al. 200li 

Iaria et alj 2004:; 

Di Salvo et al. 

2006; Paizis et al. 2006; D’Af et al. 

2007; Tarana et alJ 2007; Piraino et al. 2007). However the proper- 


ties of the hard tails are different according to the spectral state in 
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which they are observed. In particular, during the low/hard states, 
they dominate the broadband spectrum and peak above 100 keV, 
showing a cut-off at high energy which suggests an underlying 
thermal process (e.g. lMcClintock & Remillard l2006h . On the other 
hand, during the soft states, the hard tails are weak and carry only 
a few percent of the total flux emission and they do not present 
high energy cut-off. These can be explained as either being pro¬ 
duced in a hybrid thermal/non thermal medium dPoutanen & Coppil 
1 1998) or by the bulk motion of accreting material close to the NS 
(e.g. Titarchuk & Zanniatll 19981) . An alternative scenario proposes 
a mechanisms of Comptonization of seed photons by relativistic, 
non thermal, electrons in an outflow (e.g. iDi Salvo etai]|200Cl) 
or synchrot ron emi ssion from a relativistic jet moving out from 
the system jMarkoff et al.ll200lh. In the latter case, radio emissio n 
should be observed (e.g. lMigliari et al]|2007l : iHoman et alj|2004l) . 
According to the fact that GX 3+1 is more likely in a high/soft 
state and the powerlaw is weak and does not show any spectral 
cut-off, we suggest its hard tail is similar to those found during the 
high/soft state of the Galactic accreting systems._However, no radio 
emission has still been observed in GX 3+1 feerendsen et all200d) 
and a combined campaign of X-ray and radio observations would 
be required to clarify the accretion mechanisms. 

In addition to the continuum components, GX 3+1 is charac¬ 
terized by four broad emission lines which are associated to rela¬ 
tivistic reflection of hard photons from the surface of the disc. We 
inferred an inclination angle of system and a logarithm of the ion¬ 
ization parameter of the accretion disc consistent with 35° and 2.8, 
respectively. The inferred ionization parameter is in agreement with 
the H-like nature of the discrete features of the S, Ar, Ca and Fe 
that we found in the GX 3+1 spectrum. The accretion disc is more 
likely responsible for the reflection component and the shapes of 
the emission lines are strongly dependent on the underlying contin¬ 
uum, therefore a good modelling of it can support more robust con¬ 
clusions about the spectral interpretation. We argue that the study 
of the relativistic emission features of GX 3+1 with simple phe¬ 
nomenological models (e.g., DISKLINE) provides misleading inter¬ 
pretations of the accretion processes. In fact we found that the inner 
disc radius was consistent with ~ 30 R g , suggesting that the disc 
is truncated away from the NS: however this is difficult to recon¬ 
cile with a disc temperature of ~ 0.5 — 0.6 keV, as well as with 
the lack of X-ray pulsations probably reflecting the low magnetic 
field. The self-consistent models (both RFXCONV and REFLIONX) 
instead point towards more physical results: the relativistic reflec¬ 
tion is produced at a radius of ~ 10 R g which is now highly con¬ 
sistent with the disc temperature of 0.5-0.6 keV (or 0.8 — 0.9 keV, 
applying a color correction factor of 1.7). Furthermore, the disc and 
the reflection components contribute, together, to 10% of the total 
flux. We also finally note that the Ar and Ca lines in GX 3+1 are 
possibly over-abundant with respect to the solar one and further 
analyses and investigations are needed to better understand their 
nature. 

We therefore conclude that the broadband spectral properties 
of GX 3+1 are complex and well described by a combination of a 
hot accretion disc and optically thick and thin coronae likely the 
sources of the hard photons responsible for the strong reflection 
component. 

Although robust, our analysis was limited by the fact that 
the XMM-Newton and INTEGRAL datasets were not simultaneous. 
Therefore we strongly encourage further investigation of GX 3+1 
with the support of broadband X-ray satellites such as Suzaku and 
NuSTAR in order to provide a confirmation of the interpretations of 
the spectral results. 
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